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Abstract. We present new methods for determining polynomials in the ideal of the variety 
of bilinear maps of border rank at most r. We apply these methods to several cases including 
the case r = 6 in the space of bilinear maps C 4 x C 4 ->• C 4 . This space of bilinear maps 
includes the matrix multiplication operator M2 for two by two matrices. We show these newly 
obtained polynomials do not vanish on the matrix multiplication operator M2, which gives a 
new proof that the border rank of the multiplication of 2 x 2 matrices is seven. Other examples 
are considered along with an explanation of how to implement the methods. 



o 

CN 



CO 

I Acknowledgements. We thank Peter Biirgisser for important discussions and suggestions. 

^ ■ 

Lower bounds in complexity theory are considered difficult to obtain. We describe a new 
method for obtaining lower bounds on the border rank which is based on a new way to find 
polynomials that vanish on bilinear maps T : C a x C b — > C c of low border rank. 

0\ ; 1.1. Rank and border rank. Let C a * := {/ : C a — > C | / is linear} denote the dual vector 

space to C a . That is, if an element of C a is represented by a column vector of height a, then 



\ C a * corresponds to row vectors, and the evaluation is just row-column matrix multiplication. A 

bilinear map T: C a x C b —¥ C c has rank one if there exist a € C a *, /3 £ C b *, and c 6 C c such 
■ that T{a,b) = a{a)(3{b)c. The rank one bilinear maps are in some sense the simplest bilinear 

CO ' maps, and T is said to have rank r if r is the minimum number of rank one bilinear maps 

which sum to T. This r is sometimes called the tensor rank of T. If one views multiplication 
by constants as a "free" operation, then the rank differs at most by a factor of two from the 
minimal number of multiplications of variables that is needed to compute T, see [3, Ch. 14] for 
5— i ■ more information. 

Since the set of all bilinear maps C a x C b -> C c is a vector space of dimension abc, it is 
natural to talk about polynomials on the space of bilinear maps C a x C b — > C c . Unfortunately, 
one cannot test directly for the tensor rank by the vanishing of polynomials, since the common 
zero locus of the set of all polynomials vanishing on the set of bilinear maps of rank at most r is, 
typically, larger than the set of bilinear maps of rank at most r. This may be described precisely 
using the language of algebraic geometry: for the purposes of this article, we define an algebraic 
variety (or simply a variety) to be the common zero locus of a collection of polynomials that is 
irreducible, in the sense that it cannot be written as a union of two zero loci. A (proper) Zariski 
closed subset of a variety X is the common zero locus of a collection of polynomials restricted 
to X, and a Zariski open subset is the complement of a Zariski closed set. The set of bilinear 
maps of rank r is a Zariski open subset of the algebraic variety formed from the zero set of all 
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the polynomials having the set of bilinear maps of rank at most r in their zero set. The border 
rank of a tensor T is defined to be the smallest r such that all polynomials vanishing on the set 
of bilinear maps of rank at most r also vanish at T, and one writes R(T) = r. In this case, T is 
arbitrarily close, in any measure, to a bilinear map of rank r (including the possibility that the 
rank of T is r). We let cr r . a b c denote the set of bilinear maps of border rank at most r. It is an 
algebraic variety. When a, b,c are understood from the context, we simply write o~ r . 

1.2. Results. We introduce a new technique based on machine learning that finds, with high 
probability, where equations which vanish on the variety of bilinear maps of border rank at 
most r can be found. Once one knows where to look, we can use methods which began in [16] 
and were refined in [1, 5] to find the actual equations. The vanishing of these equations can be 
sometimes proved rigorously and other times proven with extremely high probability depending 
on the application. With these equations and given a tensor T, one can then show its border 
rank is greater than r if such a polynomial P satisfies P{T) ^ 0. Of special interest in this paper 
will be border rank of the matrix multiplication tensor 

2 

M 2 := e M ® e h k ® efc >* £ C 4 <g> C 4 £x) C 4 , 

i ,j,k=l 

where (e^y) is the standard basis of C 2x2 = C 4 . 

The technique we present can be applied to other implicitization problems. That is, we want 
to consider a variety 

(1) X:=W) 

where Y is (possibly a Zariski open set of) a variety and g is a system of rational functions defined 
on Y. In the bilinear case, the tensors of rank at most r is a dense subset of the algebraic set 
of tensors of border rank at most r where each tensor of rank at most r can be written as 
a sum of r tensors of rank at most one. In this case, Y is simply the product of r copies of 
the variety of tensors of rank at most one and g is the linear map corresponding to taking the 
sum. Another specific application arising in physics is the analysis of vacuum moduli space in 
(supersymmetric) field theories [8] which arise as the closure of the image under a polynomial 
map of an algebraic set. 

Before summarizing the main results, we note that a result "with extremely high probability" 
means the result was proven using numerical techniques. When we write "equations," we mean 
that we have written down the explicit equation. The main results presented here using our new 
technique are as follows. 

• No polynomial of degree less than 19 vanishes on a"6;4,4,4- 

• The degree of o"6;4,4,4) with extremely high probability, is 15,456. 

• There is a 64 dimensional space of degree 19 equations that vanish on 4 with ex- 
tremely high probability. One such equation does not vanish on the matrix multiplication 
tensor M 2 , giving a new, short proof with extremely high probability that the border 
rank of M 2 is seven. 

• We construct an equation of degree 20 that vanishes on 0"6;4,4,4> This equation does not 
vanish on M 2 , giving a new, short proof that the border rank of M 2 is seven. 

• The degree of the codimension six variety 0"6;3,4,6; with extremely high probability, is 
206,472 and its ideal is empty in degree 14. 

• The degree of the codimension two variety <J7 ; 4 4 5, with extremely high probability, is 
44,000 and its ideal is empty in degree 56. 

• The degree of the hypersurface <tb;3,5,7 is 105, with extremely high probability. 
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• The degree of the hypersurface <Ti8 ; 7,7,7 ; with extremely high probability, is at least 187,000. 

1.3. Other methods for finding equations. Very little is known about the equations of a r 
in general. One can reduce to the case ofa = b = c = r via a process called inheritance. 
Additionally, there is a systematic way to determine the equations in any given degree using 
multi-prolongation. For a discussion on inheritance and multi-prolongation, see [15, §3.7]. Even 
though multi-prolongation is systematic, it is very difficult to utilize except in very small cases. 
Most known equations have been found by reducing multi-linear algebra to linear algebra. See 
[13, 17] for the most recent equations that go up to 0"2m-2;m,m,m- Some information about the 
ideal of be found using representation theory (via the algebraic Peter- Weyl Theorem) 
as this case is an orbit closure, see [4] for an exposition. By inheritance one could deduce the 
o~r;m,n,p case for any m,n,p from the er r;T . jT . jT . case. 

1.4. Polynomials on vector spaces. We write I(cr r ) for the set of all polynomials vanishing 
on ci r , which forms an ideal. Since a r is invariant under re-scaling, we may restrict our attention 
to homogeneous polynomials since, in this case, a polynomial will be in the ideal if and only if 
all of its homogeneous components are in the ideal. 

Let V be a vector space. A subset X C V is called an algebraic set if it is the common zero 
locus of a collection of polynomials on V. Recall that we say that an irreducible algebraic set 
is a variety. If X C V is a variety that is invariant under re-scaling, let S d V* be the space of 
homogeneous polynomials of degree d on V and Id(X) C S d V* be the component of the ideal 
of X in degree d. 

Roughly speaking (see §2 for more details), our technique for studying the equations that 
vanish on a variety X of positive dimension is by applying numerical algebraic geometry and 
machine learning techniques to finite subsets of X which lie in a common linear space. That is, 
we aim to study finite subsets of algebraic sets of the form Y = X n £ C C where C is a general 
linear space of codimension at most dim X. If the codimension of C is dim X, then Y consists 
of deg X points. If the codimension of £ is strictly less than dim X, then Y is also a variety 
with the same degree as X. Moreover, if one considers X C V and Y C C, and defines dx and 
dy to be the minimal degree of the nonzero polynomials in I{X) and I(Y), respectively, then 
dx > dy- In particular, dim Id(Y) > dim Id(X) for any d < dx with similar bounds for all 
d > that can be developed from the corresponding Hilbert functions. Once we have inferred 
information about polynomials in I(Y), we use representation theory to identify which modules 
could appear. Finally, sample vectors from these modules are used to test if the entire module 
is in the ideal I{X) or not. 

2. Deciding where to go hunting 

The basic idea of our algorithm is to combine the ability of numerical algebraic geometry 
to compute points on certain subsets of a variety with methods in machine learning to obtain 
information about this subset from the computed points. We first describe needed concepts 
from numerical algebraic geometry and then from machine learning. 

At a basic level, the algorithms of numerical algebraic geometry (see [21] for general back- 
ground information) perform numerical computations on varieties where each variety is repre- 
sented by a data structure called a witness set. Let / be a polynomial system. The common 
zero locus of / is an algebraic set that can be decomposed uniquely into finitely many varieties, 
none of which is contained in the union of the others. If X is one of these varieties, called an 
irreducible component of the zero locus of /, then a witness set for X is the triple {f,L,W} 
where the zero set of L defines a general linear subspace of codimension equal to the dimension of 
X and W is the intersection of X with this linear subspace defined by L. Given one point in W, 
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arbitrarily many points on X can be computed in a process called sampling. In numerical terms, 
computing a point "on" a variety means that we have both a numerical approximation of the 
point along with an algorithm that can be used to approximate the point to arbitrary accuracy. 

This witness set description is not useful for the problems at hand since, for each of the 
varieties X under consideration, we do not assume that we have access to a polynomial system 
/ let alone any nonzero polynomials which vanish on X. However, we do assume that we have a 
description of X in the form (1). In fact, by adding variables and clearing denominators, we can 
assume that X := tt(Z) where ir is a projection map and Z is an irreducible component of the 
zero locus for some polynomial system F. This is demonstrated in the following simple example. 

Example 2.0.1. The set X := {(x, y) € C 2 | x 2 + y 2 = 1} is equal to g(Y) where 



We also have X = tt(Z) where tt(x, y, t) = (x, y) and Z is the zero locus (which is irreducible) of 



With this setup, we utilize a pseudowitness set [9, 10] for X which is the quadruple {F, ir, L, W} 
where L defines a linear subspace of codimension equal to the dimension of Z and W is the 
intersection of Z with this linear subspace defined by L. In this case, the linear polynomials L 
are constructed so that it has exactly dim X general linear polynomials in the image space of 
7r, i.e., intersect X in degX many points, while the remaining linear polynomials are general. 
In particular, n(W) consists of exactly degX distinct points. As with traditional witness sets, 
one can sample and perform membership tests on X [10]. 

The key here is that once a single sufficiently general point is known on X, other points on X 
can be computed as well. In fact, these other points can be forced to live in a fixed general linear 
subspace of codimension at most dim X thereby simplifying the future computations since one 
can work intrinsically on this linear subspace. If the intersection of the linear subspace and X is 
positive dimensional, then it is also a variety and arbitrarily many points can be sampled from 
this variety. If the intersection is zero-dimensional, it consists of exactly degX points which, 
after computing one, random monodromy loops [19] could be used to try to compute the other 
points. The trace test [20] provides a stopping criterion for deciding when exactly deg A points 
have been computed. 

Clearly, any polynomial which vanishes on X must also vanish on a finite subset of X. Al- 
though we will not delve too deep into the theory here, one can recover the invariants of X from 
a general linear subspace section of X when X is an arithmetically Cohen- Macaulay scheme (see 
[18, Chap. 1]). Nonetheless, since our current focus is on developing a list of potential places 
of where to look to focus further representation theoretic computations, we can consider all 
varieties and not just the arithmetically Cohen-Macaulay ones. Of course, this is at the expense 
of bounds rather than equality as demonstrated in the following example. 

Example 2.0.2. Consider the following varieties in P 3 : 



It is easy to verify that 

• dim X\ = 1, deg X\ = 3, and I{X\) is generated by three quadratics; 

• dim X<i = 1, deg X^ = 4, and /(A2) is generated by a quadratic and three cubics. 





X 1 :={(*V 



t, st 2 ,t 3 ) I (s, t) G P 1 } and X 2 := {(s 4 , s 3 



t,st 3 ,t 4 ) I (s,t) € P 1 }. 
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Let Yi = Xi fl % be the set of deg points where % is the hyperplane defined by the vanishing 
of £{x) = xq + 2x\ + 3^2 + 5^3. If we consider Yi C H, then 

• I(Yi) is generated by three quadratics; 

• KY2) is generated by two quadratics. 

To summarize, X\ is the twisted cubic curve in P 3 which is arithmetically Cohen-Macaulay so 
that, for example, the dimension of Id(Xi) can be determined from /^(Yi). However, X2 is not 
arithmetically Cohen-Macaulay which, in this case, can be observed since 2 = dim 12(^2) > 
dim 72(^2) = 1- Even though one should only expect dx 2 > dy 2 , we do have equality in this 
case, namely dx 2 = dy 2 = 2. 

Once we have decided on our first finite set to consider, the next task is to compute polynomials 
that vanish on this finite set, e.g., via [7]. This first computation alone provides some restrictions 
on which polynomials can be in I(X). Nevertheless, we also consider what happens when we 
add new points to our finite set. For a particular degree, there are two possible choices: either 
the originally computed polynomials will vanish at the new points or the dimension of the set 
of polynomials that vanish at all the points will decrease. In the former case, we can then move 
on to searching for higher degree polynomials not generated by these polynomials. In the latter 
case, we continue adding new points. If no polynomials of a particular degree, say d, vanish 
on some finite set, then we know that dim Id(X) = and dx > d. Thus, we try again by 
considering polynomials of degree d + 1. 

Variations of this approach can be to consider sampling points from the intersection of X 
with linear spaces of increasing dimension to see how the dimension of the vanishing polynomials 
change as less restrictions are placed on the sample points. The key in the end is to control the 
growth of the dimension of the space of polynomials under consideration since this can become 
unwieldy quickly. In particular, this method is practical for varieties X of low codimension 
since we can work implicitly on linear spaces of low dimension. Also, as the degrees of the 
polynomials under consideration increase, preconditioning becomes essential to perform reliable 
computations, e.g., see [7]. 

When the codimension is one, X is a hypersurface so that the degree of X is equal to the 
degree of the polynomial defining X. In this case, one can simply compute a pseudowitness set 
to compute its degree rather then use this interpolation based approach. For example, such an 
approach was used in [2] for computing the degree of implicitly defined hypersurfaces, which 
arise as the algebraic boundaries of Hilbert's sums of squares cones of degree 38,475 and 83,200. 

Example 2.0.3. In Example 2.0.2, we considered finite sets obtained by intersecting the curves 
with a particular hyperplane. We now use this information to limit our focus when we add 
other points to our finite set. In four variables, there is a ten-dimensional space of homogeneous 
polynomials of degree 2, but with our previously computed information, this has already been 
reduced to a seven and six dimensional space for X\ and X%, respectively. More specifically, 
the four dimensional space arising from the linear polynomial £(x) along with the three and two 
dimensional spaces, respectively, from ijjC^i) an d -^(l^)) namer y 




xo£(x),xi£(x),x\X2 + 2^1X3 + 3x2x3 + 5x§, 

X2^(x), X$£(x), x\ — X1X3, x\ — X1X3 — X2X3 — 10x| 

xq£(x), x\£(x), xix 2 + 2xix 3 + 3x 2 x 3 + 5x|, ) 

X2^(x), X$i(x), x\ — x\ + IIX1X3 + 2X2X3 + 20x| J 




By selecting additional random points, one indeed finds dim ^(-^l) = 3 and dim IiiX^) = 1. 
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This procedure develops ideas on where the degrees dj in which generators of the ideal appear. 
The next step is to conclusively determine the linear subspace of the space of polynomials of 
degrees dj that are in the ideal. For this, one uses representation theory as we describe next. 

3. Polynomials on the space of bilinear maps 

3.1. Tensors. In order to explain the polynomials it will be useful to work more invariantly, 
so instead of C a ,C b etc., we write A, B etc. for complex vector spaces of dimensions a, b etc.. 
It will also be useful to introduce the language of tensors. A bilinear map A* x B* — > C may 
also be viewed as a tri-linear map A* x B* x C* — > C, as well as in numerous other ways. To 
avoid prejudicing ourselves, we simply write T E A®B®C for any of these manifestations and 
call T a tensor. Just as we may view a linear map as a matrix after fixing bases, such T may be 
viewed as a three-dimensional matrix after fixing bases. Note that A®B®C, the set of all such 
tensors, is a vector space of dimension abc. More generally, given vector spaces A\,...,A k , 
one can define the space of tensors A±® ■ ■ ■ ® A k . There is a natural map A\® ■ ■ ■ ® A k x 
£i® ■■■® -Bj -> Ai®---®A k ®B 1 ®---®B l , (f,g) H- f®g, where f®g(an, . . . , a k , Pi, ■ ■ ■ , A) := 
f(ai,...,a k )g(p 1) ...,Pi). 

3.2. Remarks on the theory. We briefly review the representation theory underlying in the 
algorithm. For more details, see [15, Chap. 6]. Let S d (A§<)B(g>C)* denote the vector space 
of homogeneous polynomials of degree d on A®B®C. The variety <7 r;ai b,c is mapped to itself 
under changes of bases in each of the vector spaces and thus if we have one equation, we can 
obtain many more by changing bases. That is, let GL(A) denote the set of invertible linear 
maps A — »■ A and similarly for B,C. The group G := GL{A) x GL{B) x GL{C) acts on 
A®B®C by (gA,gB,gc) ■ (Yli a i®bi(&Ci) = Yli9Aai®gBbi®gcCi, and GL(V) acts on S d V* by 
g • P[x) = P(g~ 1 ■ x). Letting V = A(£>B®C and noting G C GL(V), we have a G-action on 
S d (A®B(3C)*. If P e I(cr r ), then g ■ P G I{a r ) for all g G G. Since ideals are in particular 
vector spaces, the linear span of the orbit of P in S d (A®B£>$)C)* will be in I(o~ r ). 

We will use the action of the group G to organize our calculations. A group G is said to act 
on a vector space V if there is a group homomorphism p : G —?■ GL(V). Then V is called aG- 
module. The G-module V is said to be irreducible if there is no nontrivial subspace of V invariant 
under the action of G. The irreducible polynomial GL(V) modules are indexed by partitions 
7r = (pi, . . . ,p v ), where p\ > ■ ■ ■ > p v > 0. We write |-7r| = p± + • • • + p v , and we say tt is a 
partition of |7r|. Let S n V denote the corresponding irreducible G-L(V)-module. It occurs in V r ®' 7r ' 
and no other tensor power, however not uniquely - there is a vector space's worth of realizations 
except in the cases tt = (d) or tt = (1,...,1). The irreducible GL(A) x GL(B) x GL(C)- 
modules are all of the form Va®Vb®Vc where Va is an irreducible GL(^4)-module etc.. For 
G = GL{A) x GL{B) x GL(C), every G-module decomposes into a direct sum of irreducible 
submodules. This decomposition is not unique in general, but the isotypic decomposition, where 
all isomorphic modules are grouped together, is. 

We are interested in the homogeneous polynomials of degree say d on A®B®C } denoted 
S d (A®B<g)C)* . Via polarization, a polynomial may be considered as a symmetric tensor so 
S d (A®B®C)* C {A®B®C)*® d ~ A*® d ®B*® d ®C*® d . Thus, the isomorphism types of irre- 
ducible G-modules in S d (A®B&C)* are described by triples (tt, /j,, v) of partitions of d whose 
number of parts £(n) < a etc.. Let k n ^ tU denote the multiplicity of S n A*&iS^B*®S u C* in 
S d (A<g>B®C)*, that is, the dimension of the space of realizations of S 7T A*®S tl B*®S u C* in 
S d {A®B®C)*. The inte gers k^^^i/ are called Kronecker coefficients and can be computed com- 
binatorially. The programs SCHUR or Sage or several other ones will compute them for you in 
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small cases. We used a program written by Harm Derksen, which is based on characters of the 
symmetric group. 

There is a simple formula for dim S^A* , namely 

a- a a* it Pi-Pj+j-i 

aim.D( pij ... )Pe )A — LLi<i<j<k r^-. 

see, e.g., [6, Thm 6.3]. We will be interested in cases where the dimension is small. 

Let &d denote the group of permutations of d elements. If a = b = c, then ov,a,a,a is also 
invariant under the S3-action permuting the vector spaces. Thus anytime Si Tl A*<giS n2 B*®S 1T3 C* 
is in the ideal of a r , the module S^m A*®S- Kiy{2) B*®S. Kiy ( 3) C* will be as well, for any a S ©3. 

3.3. First algorithm: to obtain a sample collection of polynomials. What follows is 
an algorithm to compute a basis of highest weight vectors for each isotypic component in 
S d (A®B®C)* . Once one has these, for each isotypic component, one can test if there are 
modules in the ideal of <r r;ai b,c (or any G-variety for that matter) by sampling random points 
on ci r;a! b,c as described in the second algorithm. For a € &d, we write cr(ui<8> ■ ■ ■ ® vj) := 
u <r(i)® "■ " ® v a(d)- Once and for all fix bases a , ...,a a of A*, and similarly for B,C. Let 
7r = (pi, . . . ,p£, 0, . . . , 0) be a partition as above. Write ir = (pi, . . . ,pi) and £(ir) = i. Define 
Fa,k G A*® d by ' 

Fa,-k ■= (a 1 )® (pi ~ P2) (8)(a 1 A a 2 )® (p2 " p % • • • (a 1 A • • • A a f )^ p f~ p f-'\ 
Here v x A • • • A v k := ^ E CTG 6 fc sga(a)a(v 1 ® ••■(g) v k ). 
Input: Degree d and partitions ir, fi, v of d. 

Output: A basis P of the highest weight space vector of the isotypic component of S n A* ®S^B* ®S U C* 
in S d (A®B®C)*. 
1: Use your favorite method to compute k-^^ v . 
2: Set k = 0. 
3: while /c < k^^^v 

do 

4: repeat 

5: Choose permutations t±,T2 € 6^. 
6: Define 

== ^®(ri • F BiAt )®(r 2 ■ F C)I/ ) G A*®*®^®*®^, 

rearrange the factors so it is expressed as an element of (A(g>B(g)C)*® d , and symmetrize 
to get 

pTl,T 2 . _ . pT!,T 2 

= X) ( ff ■ • n • F B)M )®(<7 • T 2 • ify,) G A* m ®B*®%C*® d 

<r£&d 

where recall cr • (a 1 ® • • • ® a d ) := a CT W® ■ ■ • ® a CT(d ). 
until Pn^l is linearly independent of Pi, ... , Pfc-i- 
Increase = & + 1. 
9: Set P fe = P£$. 
10: end while 



3.4. Examples. 
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3.4.1. d = 2, (ir, //, v) = ((2), (1, 1), (1, 1)). Here &(2),(i,i),(i,i) = 1? so we are looking for a single 
polynomial. We have -^4,(2) = ( fll ) 2 > Pb,(i,i) = b 1 Ab 2 and = c 1 A c 2 . Try n = r 2 = Id, 
then 

= (a 1 ®6 1 (g)c 1 )(g)(a 1 (g)6 2 (g)c 2 ) - (a 1 ®6 1 (g)c 2 )(g)(a 1 (8)6 2 (g)c 1 ) 
- (a 1 ®b 2 ®c 1 )®(a 1 ®b 1 ®c 2 ) + (a 1 ®6 2 «>c 2 )®(a 1 «)6 1 «)c 1 ) 

Thus 

tf$!{li),(il)P jka i® b 3® c k) = 2x ul x 122 - 2x n V 21 
Here, and throughout, repeated indices are to be summed over. Note that if T = x^ k ai®bj®Ck 
has rank one, then P(o)(i 1) (1 1) CO = but P w * n evaluate to be nonzero on a general rank 
two tensor. 

3.4.2. d = 3, (-7T, /i, v) = ((2, 1), (2, 1), (2, 1)). Here £^2,1), (2,1), (2,1) = so a S a i n we are looking for 
a single polynomial. We have -^4,(2,1) = a 1 <8>(a 1 A a 2 ), and similarly for B,C. Try n = T2 = Id, 
then 

F (2iU2i),(2i) = (a 1 ®^®^ 2 - a 1 ®^®^)®^ 1 ®^ 1 ®^ 2 - fe 1 ®^ 2 ®^ 1 )®^ 1 ®^®^ - c 1 ®^®^) 
= (a 1 ®b 1 ®c 1 )®(a 1 ®b 1 ®c 1 )®(a 2 ®b 2 ®c 2 ) - (a 1 (g)6 1 (g)c 1 )®(a 1 (g)6 1 (g)c 2 )(g)(a 2 (g)6 2 (g)c 1 ) 

- (a 1 (g)6 1 (g)c 1 )(g)(a 1 (g)6 2 (g)c 1 )(g)(a 2 ®6 1 (g)c 2 ) + (a 1 ®b 1 ®c 1 )®(a 1 ®b 2 ®c 2 )®(a 2 ®b 1 ®c 1 ) 

- (a 1 (g)6 1 (g)c 1 )(g)(a 2 (g)6 1 (g)c 1 )(g)(a 1 ®6 2 (g)c 2 ) + (a 1 ®b 1 ®c 1 )®(a 2 ®b 1 ®c 2 )®(a 1 ®b 2 ®c 1 ) 
+ (a 1 (g)6 1 (g)c 1 )(g)(a 2 (g)6 2 (g)c 1 )(g)(a 1 (g)6 1 (g)c 2 ) - (a 1 ®6 1 (g)c 1 )(g)(a 2 (g)6 2 (g)c 2 )(g)(a 1 (g)6 1 (g)c 1 ) 

Thus P(2\)(2i) (21) = so we need to try different n, r 2 . Take n = Id and T2 = (12). Then 
F (2iM2i),(2i) = (a l m l m 2 - a 1 ®a 2 ®a 1 )®(b 1 ®b 1 ®b 2 - b l ®b 2 ®b l )®{c l ®c l ®c 2 - c 2 ®c 1 ®c 1 ) 

= (a 1 ®b 1 ®c 1 )®(a 1 ®b 1 ®c 1 )®(a 2 ®b 2 ®c 2 ) - (a 1 ®b 1 ®c 2 )®(a 1 ®b 1 ®c 1 )®(a 2 ®b 2 ®c 1 ) 

- (a 1 ®b 1 ®c 1 )®(a 1 ®b 2 ®c 1 )®(a 2 ®b 1 ®c 2 ) + (a 1 ®b 1 ®c 2 )®(a 1 ®b 2 ®c 1 )®(a 2 ®b 1 ®c 1 ) 

- (a 1 ®b 1 ®c 1 )®(a 2 ®b 1 ®c 1 )®(a 1 ®b 2 ®c 2 ) + (a 1 ®b 1 ®c 2 )®(a 2 ®b 1 ®c 1 )®(a 1 ®b 2 ®c 1 ) 
+ (a 1 ®b 1 ®c 1 )®(a 2 ®b 2 ®c 1 )®(a 1 ®b 1 ®c 2 ) - (a 1 ®b 1 ®c 2 )®(a 2 ®b 2 ®c 1 )®(a 1 ®b 1 ®c 1 ). 

Thus i^ ( 5 ))(21) J] x t3k ai®bj®c k 

111 111 222 1 r, 112 121 211 111 121 212 111 211 122 111 112 221 

X _l_ t rp rp rp rp rp rp rp rp rp rp rp rp 

Note that if T has rank one, then f^i) (21) (21) ("0 = ^' ^ wm evama t e to be nonzero on 
a general rank two tensor. 

3.4.3. Permutation pairs to avoid. We want to avoid the case that occurs in Example 3.4.2, 
i.e., that Pn^u = 0. Although a complete classification of the cases when this happens is 
unknown, an easy necessary condition for P^ffi 7^ is the following [11, Lemma 7.2.7]: When 
we write 1,2,... ,d in a tableau columnwise starting with the longest column and we write 
ti(1),ti(2), . . . ,T±(d) in a second tableau columnwise, and we do the same for T2 in a third 
tableau, then it is necessary that there exists no pair of numbers that lies in the same column 
in all three tableaux. If this occurs, we call this situation a zero pattern. We can choose random 
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permutations that avoid the zero pattern by just choosing random permutations and repicking 
if it contains a zero pattern. 

3.4.4. d = 20, (7r,fi,i/) = ((5555), (5555), (5555)). Here £(5555), (5555), (5555) = 4. The following 
choices of pairs t±,T2 give 4 linearly independent polynomials. 

Ti = (n(l), n(2), . . . , n (20)) = (10, 15, 5, 9, 13, 4, 17, 14, 7, 20, 19, 11, 2, 12, 8, 3, 16, 18, 6, 1), 
r 2 = (10, 11, 6, 2, 8, 9, 4, 20, 15, 16, 13, 18, 14, 19, 7, 5, 17, 3, 12, 1) 

n = (19, 10, 1, 5, 7, 12, 2, 13, 16, 6, 18, 9, 11, 20, 3, 17, 14, 8, 15, 4), 
r 2 = (10, 5, 13, 6, 3, 16, 11, 1,4, 18, 15, 17, 9, 2, 8, 12, 19, 7, 14, 20) 

n = (16, 20, 9, 13, 8, 1,4, 19, 11, 17, 7, 2, 14, 3, 6, 5, 12, 15, 18, 10), 
r 2 = (1, 20, 11, 19, 5, 16, 17, 2, 18, 13, 7, 12, 14, 10, 8, 15, 6, 9, 3, 4) 

n = (11, 5, 2, 1, 16, 10, 20, 3, 17, 19, 12, 18, 13, 9, 14, 4, 8, 6, 15, 7), 
r 2 = (1, 6, 15, 13, 20, 3, 18, 11, 14, 2, 9, 5, 4, 17, 12, 8, 19, 16, 7, 10) 

This gives rise to 4 polynomials /i, • • • , /4, and the following linear combination vanishes on a§: 
-266054 /i + 421593 / 2 + 755438 f 3 + 374660 / 4 . 

3.4.5. d = 19, (-7T, /j,, v) = ((5554), (5554), (5554)). Here ^(5554)^5554)^5554) = 31. We found 31 
pairs Ti,r 2 that result in 31 linearly independent polynomials by choosing T\ and r 2 randomly, 
but avoiding the zero pattern. 

3.5. Second algorithm: to test on the secant variety. Once one has a basis of highest 
weight vectors for an isotypic component, one needs to determine which linear combinations of 
basis vectors vanish on a r . The following algorithm is standard linear algebra: 

Input: The output of first algorithm for some (tt, ji, is), i.e., a collection P\, . . . = P^ v £ 
S d (A<g>B®C)* and r, where we will test for polynomials in /(oy.at, c)- 

Output: with probability as high as you like the component of /(o" r;a] b, c ) i n S 7T A*®S^ l B*^S u C* . 
If the component is zero, then the answer is guaranteed correct, and more generally, the 
algorithm can only overestimate the component if the points on a r are not chosen randomly 
enough. 

1: Set P = c\P\ + • • • CfcPfc, where c\, . . . , c& are variables. 
2: Chose "random" vectors 

a b c 

V J = J2J2J2( a i^®(Pij b ^®^ l i,jCi) + ■■■ + (4 )i a i )(8)(/3^-6 fe )(8)(7^c z ) 
i=i k=i 1=1 

where the a Sj'^5j^8j are "random" numbers. 
3: Evaluate P at these k points. 

4: if there exist a solution c\, . . . , c\~ such that all the evaluations are zero then 

5: If there is a m-dimensional solution space, then with reasonable probability one has m 

copies of the module in the ideal. 
6: else 

7: No module in this isotypic component is in J((j r;a ^ c ). 

8: end if 

Remark 3.5.1. One can parametrize a T and in some cases use the explicit parametrization to 
get an exact answer, which is what we did for (ir,fi,v) = ((5555), (5555), (5555)). 
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4. Review of the original proof that M 2 c 6;4i4j4 

The essence of the proof that the rank of M2 is not six in [14] is as follows: there is a 
now standard argument due to Baur for proving lower bounds for rank by splitting a putative 
computation into two parts using the algebra structure on the space of matrices. The argument 
in [14] was to apply the same type of argument to each component of the variety consisting 
of subvarieties where the rank is greater than the border rank. The article [14] contained a 
gap in the proof that was filled in [12] but not published in JAMS because the editor was 
concerned the erratum was almost as long as the original article and the author did not see a 
way to shorten it. The gap in [14] was caused by overlooking the possibility of certain types of 
components, where the limiting 6-planes are not formed by points coming together but by some 
other unusual configuration of points. All such components of <T6;4,4,4 are n °t known explicitly, 
but the correction only used qualitative aspects of how the limiting 6-plane arose. There were 
3 basic cases, where any subset of 5 of the limit points are linearly independent, where there is 
a subset of 5 that are not, but any subset of four are, and where there is a subset of 4 that are 
not, but any subset of 3 are. In each of these cases, one is forced to have a limit taking place 
among rank one tensors in a much smaller space, which was what made the analysis tractable. 
The computations performed above provide an explicit polynomial vanishing on <76 ;4)4i4 which 
does not vanish at M2, providing a significantly shorter proof of this fact. 
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